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Q . During the first part of this review we will focus on the thermodynamics of SU (N) 

CSJ I gauge theories at finite temperature. We will present results from a calculation 

CN . of electric and magnetic screening masses for the gluons, discuss calculations of the 

critical temperature in units of the string tension and bulk thermodynamic quantities 
like the energy density and pressure. In particular, the latter calculations have now 
Q>^ I reached a stage where O(a^) cut-off effects can be controlled systematically and 

Csl ' an extrapolation to the continuum limit can be performed. In the second part we 

O ' discuss the critical behaviour of QCD with light quarks. We analyze the chiral 

transition in 2-fiavour QCD and present results on the temperature dependence of 
; hadron properties. 

I. Introduction 

r-| , Lattice studies of the thermodynamics of QCD have been performed since more than 

ten years The goal of these investigations is to understand quantitatively the 
^ ; behaviour of strongly interacting matter at high temperatures. In order to achieve 

' this computer simulations of lattice regularized QCD have to be performed with pa- 

rameters for the quark masses and number of flavours which are as close as possible 
to the values in the real world, ie. a world where the dynamics is controlled by two 
light, nearly massless quark flavours. In addition these calculations should be per- 
formed as close as possible to the continuum limit in order to minimize systematic 
errors resulting from the finite lattice cut-off. Numerical calculations are steadily 
improving in both respects. Nonetheless, given the presently available computa- 
tional resources, one still has to compromise either by performing calculations in 
the physically interesting regime with light quarks but relatively large lattice spac- 
ings and limited statistical accuracy or by analyzing the approach to the continuum 
limit through high statistics calculations within the context of quenched QCD, ie. 
in pure SU{N) gauge theories. 

Most of the results with high statistical accuracy have been obtained for pure 
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SU {N) gauge theories, which do aheady provide a highly non-trivial model for the 
thermodynamics of QCD. Many of the global features of the thermodynamics like, 
for instance, the rapid rise of the energy density in the plasma phase towards the high 
temperature ideal gas limit, have been found to be qualitatively similar in SU{N) 
gauge theories and QCD with light quarks. Moreover, it is well known that also the 
finite cut-off effects, which are quite large in calculations of bulk thermodynamic 
observables, are similar in QCD and SU{N) gauge theories. The simpler SU{N) 
gauge theory thus is an ideal model for the systematic investigation of discretization 
errors and their controlled removal in the continuum limit. We will discuss this in 
section II of these lectures. 

The investigations of the QCD thermodynamics in the presence of light quarks, 
on the other hand, are not yet on the same quantitative level. The analysis of the 
equation of state, for instance, still suffers from large discretization errors. Another 
important issue in the analysis of the thermodynamics with two light quarks is to 
understand the role of chiral symmetry restoration at the critical temperature. As 
this is an infrared problem, the discretization errors do not play a central role in this 
case. Here the main objectives are to model the correct continuum flavour symmetry 
in lattice regularized QCD and to perform calculations at small enough quark masses 
in order to be able to observe the influence of chiral symmetry breaking/ restoration 
on the QCD thermodynamics. We will discuss this aspect in section III of these 
lectures. 

In this short presentation we will not be able to review systematically the for- 
mulation of lattice regularized QCD. For a more detailed discussion of the basic 
concepts with emphasis put on the finite temperature formulation we refer to recent 
reviews and references given therein 0, 0, ^ . 

II Equation of State for SU(N) Gauge Theories 

II. 1 Screening masses and perturbation theory 

Our intuitive understanding of the high temperature phase of QCD and as such 
also of the non-abelian SU (N) gauge theories has been guided to a large extent 
by perturbative calculations at high temperatures. Asymptotic freedom suggests 
that at very high temperatures the QCD plasma approaches an ideal gas of non- 
interacting quarks and gluons. This requires that the temperature dependent run- 
ning coupling, g{T), becomes small. At the same time, it is well-known that per- 
turbative calculations are faced with severe infrared problems which only can 
be overcome, if thermal electric and magnetic masses of order g{T)T and g^{T)T, 
respectively, are generated for the gluons. Even then the perturbative expansion of 
the thermodynamic potential is, however, known to show quite poor convergence 



2 



properties also at temperatures which are several times as large as the critical tem- 
perature [0. One thus may wonder to what extent the concept of magnetic and elec- 
tric screening masses, in particular the arguments given for their relative ordering, 
me(T) ~ 0{g{T)T) > mm{T) ^ 0{g'^{T)T), are of relevance beyond perturbation 
theory. After all, such arguments rely on the assumption that g{T) will be small at 
high temperature. 

A recent analysis of the gluon propagator in a SU{2) gauge theory [0, performed 
in Landau gauge, shows that the expected ordering of magnetic and electric screening 
masses indeed persists down to temperatures close to Tc, ie. me(T) > rrimiT). 
However, it also shows that at least the leading order perturbative relations for 
the electric mass do not hold for temperatures at least up to 15Tc. Rather than 
decreasing with temperature me(T)/T does stay constant and is quite large, i.e. 
me(T)/T ~ 2.5. Although the magnetic mass does seem to scale like ~ g'^{T)T 
and is smaller than rrie above T^, also in this case the running coupling extracted 
from such a scaling ansatz turns out to be large, i.e. g'^{2Tc) ~ 4. Results from a 
calculation in Landau gauge [0 are shown in Figure [|. 
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Figure 1: Electric and magnetic gluon screening masses for the SU{2) gauge theory 
calculated in Landau gauge [5]. For the electric mass we show a straight line fit, rrie = 
(2.484 ± 0.052)T. For the magnetic mass we include the observed (7^-dependence, 
m„ = (0.466 ± 0M5)g^{T)T, with g-^{T) = (ll/127r2) ln((3.8 ± 0.3)T/T^) 



Although the gauge dependence of the results presented above still has to be an- 
alyzed carefully, they do suggest that the running coupling constant is larger than 
unity for temperatures up to a few times Tc. This observation is also in agreement 
with conclusions drawn from the temperature dependence of other observables like, 
for instance, the spatial string tension |^. One thus cannot expect that pertur- 
bation theory provides a quantitative description of bulk thermodynamics in this 
temperature regime A non-perturbative approach is needed. 
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II. 2 Cut-off dependence of thermodynamics 



Lattice calculations of the energy density (e), pressure (p) and other thermody- 
namic variables led to some understanding of the temperature dependence of these 
quantities in the QCD plasma phase. The energy density, for instance, has been 
found to rise rapidly at Tc and approach the high temperature ideal gas limit from 
below. Such calculations are performed on finite Euclidean space-time lattices of 
size A^^ X Nt- with non-zero lattice spacing a, which is controlled through the di- 
mensionless coupling j3 = 2N/g^. This fixes also the thermodynamic parameters, 
temperature T = and volume V = {N„a)^. In general such calculations suf- 

fer from two different types of finite size errors. An infra-red cut-off is introduced 
through the finite spatial size and an ultra-violet cut-off, at fixed temperature, 
is introduced through the temporal extent A^,-. While N^^ should be large in order to 
reach the thermodynamic limit {V — > oo), A^^ should be large in order to eliminate 
the influence of the regularization (finite lattice cut-off) on physical observables in 
the continuum limit. In the following we mainly will be concerned with the latter 
aspect. 

As a consequence of the O(a^) discretization errors in the standard Wilson formu- 
lation of SU{N) gauge theories thermodynamic quantities calculated on the lattice 
receive 0((aT)^) corrections. For an ideal gluon gas these corrections are known to 
be as large as 50% on lattices with A^,- = 4, ie. aT = N~^. Asymptotically, in the 
limit c/^ — i> 0, the cut-off dependence of energy density and pressure are given by 
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Cut-off effects are particularly large in observables, which are sensitive to the 
ultra-violet behaviour of the theory. This is the case for bulk thermodynamic quan- 
tities in the high temperature limit. As will be discussed in the following these 
finite cut-off effects are, in fact, clearly visible in the calculations of the pressure 
and energy density performed on lattices with varying temporal extent, although in 
the temperature regime up to a few times Tc they turn out to be about a factor 2 
smaller in magnitude than calculated analytically in the T — ^ oo limit Only 



recently these calculations could be extended to lattices with sufficiently large tem- 
poral extent (A^r = 6 and 8) that would allow an extrapolation of lattice results for 
bulk thermodynamic quantities to the continuum limit |jlO|, |12| . 

On the other hand, cut-off effects are of less importance for the analysis of 
long-distance properties of the theory, such as the string tension and the decon- 
finement transition temperature itself. In this case it is more important to reach 
the large volume (thermodynamic) limit. The critical temperature is calculated at 
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finite lattice cut-off, ie. on lattices with fixed temporal extend N^-, by determining 
the corresponding critical couplings [3c{Nt) |[T^, [T^, |T^, [T^. A calculation of the 
string-tension at this value of the coupling on a large zero temperature lattice is 
then used to set the scale, 

^) = (^) + c,{aa^) ■ (2) 

V / |a=0 V V ^ ^ \aT=l/Nr 

A collection of results for this ratio is shown for the SU{2) and SU{3) gauge 
theories in Figure |^. As can be seen the O(a^) corrections are quite small. An 
extrapolation to the continuum limit yields 

1^ ^ r 0.69 ±0.01 , SU(2) 

a l0.629±0.003(+0.004) , SU(3) ' 



where the number in brackets, given in the case of SU{3), indicates the systematic 
error still present in the data due to the missing extrapolation to the thermodynamic 
limit on the larger lattices. 
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Figure 2: The critical temperature in units of the square root of the string tension 
for the SU{2) and SU (3) gauge theories versus the lattice cut-off squared. 



Let us now turn to a discussion of bulk thermodynamic observables. The pressure 
can be obtained from an integration of the difference of action densities at zero (Sq) 
and finite (St) temperature [0, 
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where the action densities themselves may be chosen appropriately in order to reduce 
cut-off effects. The simplest choice corresponds to the standard Wilson action, which 
leads to 0((aT)^) errors in thermodynamic observables. It is given in terms of 
plaquette expectation values, 

^St? = 6(1 - ^W(^^^))o(T) , (5) 

with denoting the product of the four link- variables around an elementary 

plaquette in the lattice. Various improved actions have been suggested, which reduce 
the cut-off effects to 0(((aT)''). One possibility is to add to the single plaquette 
term in the Wilson action a contribution of (2 x 2)-loops, W^'^^'^\ (for further details 
see m). 



q((2x2) 



(6) 



Using either of these actions in the simulation one can perform the numerical integra- 
tion defined in Eq. ^. This yields the pressure difference between two temperatures, 
corresponding to the two couplings i3q and (3. The lower temperature, corresponding 
to Po, can be chosen small enough so that the pressure can be approximated by zero 
at this point. As the low temperature phase in a SU{N) gauge theory is a glueball 
gas, ie. a system with rather heavy constituents, the pressure drops exponentially 
for T < Tc. The integration therefore can be stopped rather close to T^. Results 
obtained for the pressure with the standard Wilson action on lattices with temporal 
extent A',- = 4, 6 and 8 as well as with an improved action on a Nr = 4 lattice 
are shown in Figure |^. One clearly sees the expected cut-off dependence. As dis- 
cussed above, it qualitatively reflects the A",-- dependence of the free gluon gas, which 
also is shown by horizontal lines in this figure. We note that the improved action 
calculation, performed on a lattice with only four sides in the temporal direction in- 
deed leads to a smaller cut-off dependence than the corresponding standard Wilson 
action. 

Making use of basic thermodynamic relations we can then evaluate the energy 
density in the thermodynamic limit from the trace anomaly, 

e-3p , /^4. r,i.Ti( . ^9' 



where the derivative adg~'^/da is obtained from an explicit parameterization of the 
dependence of the cut-off, a, on the bare coupling g"^ This is shown in Figure ^. 

The much slower rise of the pressure relative to the energy density results in a 
large peak in the trace anomaly, (e — 3p)/T^, at Tpeak — l-lT'c. The trace anomaly 
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Figure 3: The pressure versus T/Tc for A^.^ = 4, 6 and 8 obtained from an integration 
of the action density for the standard one-plaquette Wilson action. For A^.^ = 4 we 
also show the result obtained with an improved action (4(2x2))- 

defines, the difference between the gluon condensate at zero and finite temperature 
e-3p = G^{0) - G^{T). From Figure | we find at Tp^ak, 

(e - 3p)peak = (0.57 ± 0.02)(t2 ~ 2.3 GeV/fm^ (8) 



where we have used the result given in Eq. ^ for the critical temperature in units 
of the string tension and y/a ~ 420 MeV in order to set the physical scale. This 
result is, in fact, compatible with conventional values for the zero temperature gluon 
condensate, G^(0) ~ 2 GeV/fm^ It, therefore, suggests that G'^(T) ~ at Tpeak- 



1 1. 3 Extrapolation to the continuum limit 

Based on the analysis of the pressure and energy density on various size lattices 
one can now attempt to extrapolate these quantities to the continuum limit. At 
present this can only be done for the Wilson action where results do exist on lattices 
with temporal extent up to A^,- = 8. Using an ansatz motivated by Eq. |l|, 

p \ ( p \ ^ const. 



we extrapolate the A^.^ = 6 and 8 data, respectively [|T^. The resulting continuum 
predictions for the pressure, energy density and entropy density are shown in Fig- 
ure |[ We generally find that the difference between the extrapolated values and the 
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Figure 4: The difference (e — 3p)/T'^ obtained from calculations with the standard 
Wilson action on lattices with temporal extent A^,- = 4, 6 and 8. 

results for = 8 is less than 4%, which should be compared with the correspond- 
ing result for the free gas, where the difference is still about 8%. This suggests that 
relative to the ideal gas case more low momentum modes, which are less sensitive 
to finite cut-off effects, contribute to thermodynamic quantities. 

Knowing the energy density and pressure as a function of temperature we also 
can evaluate other thermodynamic quantities. A quantity of particular interest for 
the description of the hydrodynamic expansion of a quark-gluon plasma created in 
a heavy ion collision is the velocity of sound, 

,10, 

Already from Figure ^ it is clear that v"^ will be small close to because the pressure 
changes only little there, while the energy density rises rapidly. This, indeed, is seen 
in Figure ^. We also note that the velocity of sound rapidly approaches the ideal 
gas limit, vf^^^^ = 1/3. 

The above discussion of energy density and pressure in a SU{N) gauge theory has 
shown, that both quantities approach the ideal gas limit from below. We note that 
the energy density rapidly rises to about 85% of the ideal gas value at 2Tc and then 
shows a rather slow increase, which is consistent with a logarithmic increase as one 
would expect if the coupling g{T) changes only logarithmically with temperature. 
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Figure 5: Extrapolation to the continuum limit for the energy density, entropy 
density and pressure versus T /T^.. The dashed horizontal line shows the ideal gas 
limit. The hatched vertical band indicates the size of the discontinuity in e/T^ 
(latent heat) at T^. [10]. Typical error bars are shown for all curves. 

We have not discussed here the status of the equation of state for QCD with light 
quarks. A detailed investigation for two flavour QCD has been performed recently on 
lattices with temporal extent A^^ = 4 |T^. However, also here an analysis on larger 
lattices is still needed in order to control the obviously present large discretization 
errors. 

Ill QCD with Light Quarks 
III.l The chiral phase transition 

A well established result from lattice calculations is the strong dependence of the 
order of the phase transition on the QCD symmetries related to the colour as well 
as the flavour degrees of freedom. In the absence of quarks, ie. in pure SU{N) 
gauge theories, the phase transition is second order for N = 2 and first order for 
A^ = 3. In the case of QCD with nf light quark flavours the transition is first order 
for Uf > 3 and does seem to be continuous for Uf = 2. Lattice calculations also 
show a strong dependence of the transition temperature on the number of partonic 
degrees of freedom. The transition temperature, expressed for instance in units of 
the square root of the string tension {Tc/^/a), is substantially smaller in QCD with 
light quarks than in a pure gauge theory. In the case of QCD with two light quarks 
one finds Tc{nf = 2) ~ ISOMeV, while in the purely gluonic theory, nf = 0, Tc 
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Figure 6: Velocity of sound for the SU (3) gauge theory calculated on a lattice with 
temporal extent N^- = 6. We do not observe any significant cut-off dependence for 
this derivative. 

is found to be substantially higher, Tc(nj = 0) ~ 260MeV. As discussed in the 
previous section this may be understood in terms of a percolation threshold: In the 
absence of quarks there are no light hadrons, the lightest states are glueballs with a 
mass of 0(1 GeV). Therefore, a rather high temperature is needed to excite enough 
of these heavy glueball states in order to reach a critical density where hadrons start 
overlapping. 

With increasing temperature non-perturbative properties of the QCD vacuum, 
like confinement and chiral symmetry breaking, get modified through the thermal 
heat bath and eventually disappear above T^. No n- vanishing condensates play an 
important role in understanding these non-perturbative features of QCD. They vary 
with temperature in the vicinity of the QCD phase transition. The restoration of 
chiral symmetry at will, for instance, be reflected in a characteristic temperature 
dependence of the chiral condensate, (ipip) ~ \T — Tc\^, which also should have con- 
sequences for the properties of hadrons close to T^. Effective theories, deduced from 
QCD in order to describe low energy properties of hadrons, establish a close link be- 
tween hadronic properties and the non-perturbative structure of the QCD vacuum. 
This is, for instance, used in the operator product expansion (OPE) for correlation 



functions of hadronic currents which allows to relate hadron masses to conden- 
sates of the quark and gluon flelds. The temperature and density dependence of the 
latter has been discussed within the context of chiral perturbation theory ||l] 



We will discuss here some of the recent results of lattice calculations on the phase 
transition in two-flavour QCD as well as an analysis of hadronic properties close to 
the QCD phase transition. Our discussion is limited to the staggered formulation 
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of lattice regularized QCD. First attempts to study the chiral symmetry restoration 
also with Wilson fermions have been reviewed in In the next section we briefly 
discuss the critical behaviour in two-flavour QCD. Section III. 3 is devoted to a 
discussion of lattice calculations of hadron masses and decay constants at finite 
temperature within the quenched approximation. 



III. 2 The chiral transition in two-flavour QCD 



The nature of the chiral phase transition does seem to depend strongly on the number 
of light quark flavours, Uf. The special role of two-flavour QCD has been noticed 



already some time ago by Pisarski and Wilczek [Q. They have suggested that 
the dynamics of the chiral phase transition is controlled by an effective, 3-d scalar 
Lagrangian constructed in terms of the chiral order parameter fleld, a ~ ipip. This 
effective theory has a global 0(4) symmetry {SU (2) x SU (2)), which is spontaneously 
broken to SU (2). It therefore is expected that in the case of a continuous transition 
the critical behaviour is governed by critical exponents of a 3-d, 0(4) symmetric 
spin model ||20|, A remaining question in this scenario is to what extent the 



axial U{1) influences the critical behaviour in the vicinity of Tc [pi] , |2^ . An even 
more fundamental problem is whether an effective scalar theory can at all describe 
the critical behaviour of a theory, in which the scalar flelds are only constructed as 
fermionic bilinears. It recently has been argued by Kocic and Kogut that this may 
not be the case and that two-flavour QCD may instead show mean fleld behaviour 



in the vicinity of . In lattice calculations an additional problem arises: None 

of the fermion Lagrangians used for numerical calculations has the complete chiral 
symmetry of the continuum Lagrangian. In the staggered fermion formulation, for 
instance, one has a chiral U{1) rather than the full SU{2) symmetry, which only 
is restored in the continuum limit. This too may influence the critical behaviour 
in numerical studies. A detailed quantitative investigation of the critical behaviour 
thus is needed in different lattice formulations. Here we will restrict our discussion 
to the case of staggered fermions. 

The critical behaviour of QCD with light quarks is controlled by two dimen- 
sionless parameters, the reduced temperature t = (T — Tc)/Tc and the scaled quark 
mass, h = iriq/T. The latter corresponds to an external magnetic fleld in spin mod- 
els. In the case of QCD with dynamical quarks it is the quark mass rather than the 
spatial lattice volume, which is the most stringent limitation for large correlation 
lengths in the vicinity of the critical point. Rather than performing a flnite size 
scaling analysis, as it has been done in the case of the deconflnement transition in 
SU{N) gauge theories, it thus is more appropriate to consider the scaling behaviour 
in terms of the quark mass. Finite volume effects may in a flrst approximation be 
ignored. 

In the vicinity of the critical point the behaviour of bulk thermodynamic quan- 
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Figure 7: The chiral susceptibility in two- flavour QCD versus Q/g"^ for three values 
of the quark mass on a 8'^ x 4 lattice. 

titles is related to thermal {yt) and magnetic {yh) critical exponents, which charac- 
terize the scaling behaviour of the singular part of the free energy density, 



T 

fs{t,h) = -—InZs 



Here h is an arbitrary scale factor. Various scaling relations for thermodynamic 
quantities can be derived from Eq. |Tl| p4| , p5[| . For instance, one finds for the chiral 



order parameter, (■?/'■?/') , and its derivative with respect to the quark mass (the chiral 
susceptibility Xm) 



Xmit, h) 



h^/^F{z) 



6 



F{z) - ^F'iz) 



(12) 



with scaling functions F and F' that only depend on a specific combination of the 
reduced temperature and scaled quark mass, z = th~^/^^ . The critical exponents (3 
and 5 are given in terms of yt and yn as {3 = {1 — yh)/yt and 5 = ?//i/(l — yh)- A 
direct consequence of Eq. is, for instance, that Xm has a peak located at some 
fixed value of Zc- This defines a pseudo-critical point, tc{h) = Zch^^^^, at which the 
peak height of Xm increases with decreasing quark mass according to Eq. |12[ Fig. ^ 
shows the scaling behaviour of Xm for three values of the quark mass The 
rapid rise of Xm with decreasing quark mass can be used to determine the critical 
exponent 6. 



Another useful observable for the determination of the critical exponent 6 is the 
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chiral cumulant, 



hX,m 

{xx) 



(13) 



From Eq. it follows that A is directly a scaling function, ie. it only depends on 
z. Moreover, it is uniquely determined at z = 0, where it takes on the value 1/6 
even at non-vanishing values of the quark mass. The chiral cumulant thus allows a 
direct calculation of this critical exponent. In Fig. ^ we show a collection of the 
pseudo-critical couplings obtained from Monte Carlo calculations with various values 
of the quark mass and on various different lattice sizes. The curves show various 
fits for the quark mass dependence of these pseudo-critical couplings. Apparently 
it is difficult to distinguish 0(4) exponents from mean field exponents in this way. 
We note, however, that the location of the zero quark mass critical coupling is 
quite well determined. The cumulant A is shown in Fig. ^b. For the two smaller 
quark masses this quantity has been evaluated in the region of the estimated zero 
quark mass critical coupling, l3c{mq = 0) ~ 5.243(10). In this interval the cumulant 
takes on values between 0.21 and 0.26, which is compatible with the 0(4) exponent 
1/(5 = 0.205(45). 
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Figure 8: The pseudo-critical couplings (a) of the chiral transition in two-fiavour 
QCD on lattices with temporal extent = 4 and the chiral cumulant (b). Shown 
are results from calculations with different values of the quark mass on a 8^ x 4 
lattice. 



Further scaling relations can be obtained by considering the quark mass de- 
pendence of the thermal susceptibility, Xt ~ d{ipip)/dt, and the specific heat. A 
combined analysis of these allows a determination of the critical exponents yt and 
Uh- Results from a first analysis of this kind [25| are summarized in Table |l| and 
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compared with the known exponents for 3-d, 0(4) symmetric spin models and 
mean field behaviour. 



Table 1: Thermal {yt) and magnetic {yh) exponents from simulations of two-fiavour 
QCD on an 8^ X 4 lattice (LGT). Also given are mean field exponents (MF) and 
those for a 3-d, 0(4) symmetric spin model. 





LGT 


0(4) 


MF 


yt 


0.69(7) 


0.446(5) 


0.5 


Vh 


0.83(3) 


0.830(6) 


0.75 



As can be seen from the Table the exponent yt extracted from lattice calculations 
is in very good agreement with the expected 0(4) model behaviour. It corresponds 
to a value of 0.205(45) for the exponent 1/5, which is significantly below the value 
1/3 expected for a mean field critical exponent. The exponent yt, however, comes 
out too large in the present calculation. It may be expected that the determination 
of yt is more sensitive to finite volume effects than that of yh- In the limit of 
vanishing quark masses it is yt that controls the finite volume scaling behaviour of 
thermodynamic quantities. 

Clearly the analysis of critical exponents on a 8^ x 4 lattice can only be considered 
as a first exploratory step towards a detailed investigation of critical exponents in 
two-fiavour QCD. In the near future calculations on larger lattices with smaller quark 
masses will certainly lead to much better determinations of critical exponents. This 
should allow to distinguish between mean field and 0(4) critical behaviour. 

The chiral susceptibility, Xm, has a very narrow peak at the pseudo-critical point. 
This indicates that only in this region contributions from the singular part of the free 
energy, which is responsible for the occurrence of a phase transition, is dominant. 
This behaviour is consistent with the temperature dependence of (ipip) , which leads 
to deviations from the zero temperature condensate only close to Tg. We thus expect 
that also hadronic properties will show significant modifications only in a narrow 
temperature regime close to Tdrnq) = Tpeak, 



-^peak/2 -^peak 



peak 



<0.1 . (14) 



III. 3 Hadronic properties close to 

Information about hadron masses and decay constants can be extracted in lattice 
simulations from the long-distance behaviour of correlation functions of hadron op- 
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erators 



Gh{x) = {H{x)H+{0)) e-'"^!"' , x = {t,x) , (15) 

where H{x) denotes an operator with the appropriate quantum numbers of the 
hadronic state under consideration, for instance H{x) = ^u{,x)'~^;^ipd{,x) for the pion. 
At zero temperature one studies the behaviour of the correlation function at large 
Euclidean times r. From the exponential decay of Gh one deduces the hadron mass 
of the lightest state in this channel, whereas the amplitude is related to the decay 
constant of this hadronic state. At finite temperature the Euclidean extent is limited, 
< r < 1/T, and one thus studies the behaviour of Gh for large spatial separations, 
\x\ oo. This yields information about the finite temperature screening masses 
which are related to pole masses as long as there is a bound state in the quantum 
number channel under consideration. 

Finite temperature screening masses have been studied in lattice simulations of 



QCD for quite some time [^. A very drastic qualitative change in the screening 
masses is seen when one crosses the QCD transition temperature. Parity part- 
ners become degenerate above Tc, the pseudo-scalar mass becomes massive and 
approaches the free quark/anti-quark value, mmeson = 27rT, at large temperatures. 
These features do not seem to depend much on the number of quark flavours. In 
particular, they also have been found in quenched QCD simulations. It thus seems 
to be meaningful to first study the temperature dependence of hadronic properties 
in the quenched approximation where results on large lattices with high accuracy 
can be obtained. In the following we will describe the results of such an investi- 
gation performed on a rather large lattice (32^ x 8) close to the phase transition 
temperature |^ . 



III. 3.1 The GMOR relation and the pion decay constant 

The GMOR relation gives the chiral condensate in terms of the pion mass and pion 
decay constant, 

flml = mq{iljilj)^^=o . (16) 

Below Tc the pion is a Goldstone particle, its mass squared depends linearly on the 
quark mass, 

ml = a^rug . (17) 

A calculation of the chiral condensate and the pion mass at different values of the 
quark mass allows the determination of the pion slope, a^, and the zero quark mass 
limit of the condensate, {'ijj'ip)mq=o- 
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The pion decay constant /^r can be determined directly from the relevant matrix 
element 



V2Uml = (0|V'.75^dk+) 



(18) 



f 



,9=6. 00 
r=0 . 92T 




(a) 



f (T) /f , (T=0) 



(b) 



o Np=0 (m-O) 
A N„=4 {m=0.01) 



1.5 2.0 

T/T 



Figure 9: Test of the GMOR relation at T = 0.92 Tc in quenched QCD (a) and 
the temperature dependence of the pion decay constant (b). In (a) we compare 
the results obtained from the amplitude of the pion correlation function through 
Eq.(15) (circles) and extrapolated to zero quark mass with the result obtained from 
the GMOR relation (star, on top of the circle at m^a = 0). Also shown is the zero 
temperature result at this value of the gauge coupling (square). 

The square of the matrix element appearing on the right hand side of Eq. |18| is 
proportional to the amplitude of the pion correlation function and can be determined 
in a Monte Carlo calculation. A comparison of /^r determined this way with the ratio 
(^^)m,=o/'^7r thus provides a direct test of the GMOR relation at finite temperature. 
This is shown in Fig. ^a at T ^ 0.92Tc. Moreover, we can compare the value of 
Jtt calculated at finite temperature with corresponding zero temperature results. 
This is shown in Fig. ^d. Clearly we do not have any evidence for violations of 
the GMOR relation nor for a significant change of /^r with temperature below T^. 
The sudden change above Tc reflects the drastic change in the structure of the 
pseudo-scalar correlation function. It does no longer have a pole corresponding to 
a Goldstone-particle (pion). The existence of a pseudo-scalar bound state above Tc 
thus is questionable. Certainly for large temperatures such a state does not exist, 
the correlation function is dominated by a quark/ ant i- quark cut. The notion of f^^ 
used for the square root of the amplitude of the correlation function should thus be 
used with caution above Tc. 



III. 3. 2 Vector meson mass and nucleon mass 



16 



The possibility of a variation of the p meson mass with temperature has been dis- 
cussed a lot as this might lead to modifications of dilepton spectra, which are exper- 
imentally detectable. The temperature dependence of the meson masses has been 
discussed within the framework of the OPE. Arguments have been given that the 
meson masses are temperature independent up to O(T^) |^9|]. The Monte Carlo 
calculations of the vector meson correlation function at finite temperature also show 
no significant temperature dependence of the mass even close to T^. In Fig. |10| we 
show the result of our calculation of the screening mass in the vector channel corre- 
lation function at T 0.92Tc and compare this with zero temperature calculations 
at the same value of the gauge coupling. There is no evidence for any tempera- 
ture dependence. The same holds true for the nucleon mass, although the details 
are more subtle in this case. As can be seen in Fig. p!0[ b there is a clear difference 
between the local masses, mjq{z) \n.G n{z) / G n{z + 1) , extracted on a large zero 
temperature lattice and those extracted at finite temperature on a lattice of size 
32^ X 8. However, as the nucleon is a fermion, there is a non-negligible contribution 
from the non-zero Matsubara mode to the nucleon screening mass. 



mN = ^/m% + (7rT)2 (19) 

After removing the thermal contribution, irT, the nucleon mass, "friN, agrees with 
the zero temperature result within statistical errors. 

IV. Conclusions 

We have discussed recent detailed studies of the equation of state of the SU (3) gauge 
theory. Deviations from an ideal gas behaviour are still significant at temperatures 
of a few times Tc. The now existing extrapolation of the SU{3) equation of state 
to the continuum limit does provide a good basis for quantitative comparisons with 
perturbative calculations and phenomenological models PD|. 



The discussion of basic properties of the thermodynamics of two-flavour QCD 
shows that there are strong indications for a second order phase transition. However, 
whether two-flavour QCD is, indeed, in the same universality class as the 3d, 0(4) 
model cannot yet be decided on the basis of the existing numerical data and needs 
further detailed investigations on larger lattices and with smaller quark masses. 

At least up to temperatures of about O.QTc so far no statistically significant 
temperature dependence could be observed for basic properties of hadrons. This is 
true for the chiral condensate, hadron masses as well as decay constants calculated 
in quenched QCD. At least for the chiral condensate one expects, however, a strong 
variation with temperature for 0.9Tc < T < in particular in the case of a second 
order phase transition. This is clearly supported by the currently available lattice 
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(a) 



o g 



(b) 



,9=6. 00 
m =0. 01 



Figure 10: The p-meson mass in quenched QCD for various values of the quark 
mass (a) and local masses from the nucleon correlation function (b). The simulations 
have been performed at fixed value of the gauge couphng, 6/g'^=6.0, on a 32^x8 
lattice corresponding to T ~ 0.92 Tc and on various low temperature lattices [5]. In 
(b) we show estimates for the nucleon mass obtained from the nucleon correlation 
function at distance z. These local masses converge to the nucleon mass in the limit 
2; — > oo. The horizontal band indicates the corresponding extrapolation at zero 
temperature. No temperature dependence is visible in the p-meson mass as well 
as the nucleon mass (see text for appropriate subtraction of the finite temperature 
Matsubara mode). 



calculations in two-flavour QCD. The influence of this on the behaviour of hadronic 
parameters will certainly also be investigated in the near future. 
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